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ABSTRACT 

The direct detection of exoplanets has been the subject of intensive research in the 
recent years. Data obtained with future high-contrast imaging instruments optimized 
for giant planets direct detection arc strongly limited by speckle noise. Specific ob- 
serving strategics and data analysis methods, such as angular and spectral differential 
imaging, are required to attenuate the noise level and possibly detect the faint planet 
flux. Even though these methods are very efhcient at suppressing the speckles, the 
photometry of the faint planets is dominated by the speckle residuals. The determi- 
nation of the effective temperature and surface gravity of the detected planets from 
photometric measurements in different bands is then limited by the photometric error 
on the planet flux. In this work we investigate this photometric error and the con- 
sequences on the determination of the physical parameters of the detected planets. 
We perform detailed end-to-end simulation with the CAOS-based Software Package 
for spectro-polarimctric high-contrast cxoplanet research (SPHERE) to obtain real- 
istic data representing typical observing sequences in Y, J, H and Ks bands with a 
high contrast imager. The simulated data are used to measure the photometric accu- 
racy as a function of contrast for planets detected with angular and spectral-|- angular 
differential methods. We apply this empirical accuracy to study the characterization 
capabilities of a high-contrast differential imager. We show that the expected photo- 
metric performances will allow the detection and characterization of exoplanets down 
to the Jupiter mass at angular separations of I.O" and 0.2" respectively around high 
mass and low mass stars with 2 observations in different filter pairs. We also show that 
the determination of the planets physical parameters from photometric measurements 
in different filter pairs is essentialy limited by the error on the determination of the 
surface gravity. 

Key words: techniques: high angular resolution - techniques: photometric - methods: 
data analysis - infrared: planetary systems 



1 INTRODUCTION 

Since the detection of the first exoplanet orbiting a main 
sequence star, 51 Peg, a large population of these objects 
has been discovered covering a wide range of masses and 
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orbital periods mostly with indirect methods such as ra- 
dial velocities measurements and transits (see e.g. Santos 
2008 for a review). Although mainly sensitive to plan- 
ets with period of less than 10 years, radial velocities 
surveys have found stars which start to show long term 
trends indicating possible low mass companions orbiting 
at large orbital separations (Wittenmyer et al. 2007). The 
wide use of adaptive optics (AO) systems and coronagra- 
phy in large telescopes instrumentation for high-contrast 
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imaging has allowed to start probing the vicinity of nearby 
stars for low mass companions at large orbital separa- 
tions. Over the last decade, a handful of objects close to 
the planetary mass regime have been imaged with exist- 
ing instruments, such as 2M 1207 b (Chauvin et al. 2005), 
DH Tau B (Itoh et al. 2005), GQ Lup b (Neuhauser et al. 
2005), AB Pic b (Chauvin et al. 2005), CHXR 73 B 
(Luhman et al. 2006), and more recently Fomalhaut b 
(Kalas et al. 2008), IRXS J1609 b (Lafreniere et al. 2008), 
P Pic b (Lagrange et al. 2008) and the triple system around 
HR 8799 (Marois et al. 2008). However, the large uncer- 
tainty on the mass of these objects may place some of them 
in the sub-stellar rather than the planetary mass regime. 

The next generation of planet finding instruments cur- 
rently being built will combine: (i) extreme AO systems with 
a large number of actuators (Angel 1994; Stahl & Sandler 
1995; Langlois 2001) to reach very high corrections in 
the near-infrared (Fusco et al. 2006; AUer-Carpentier et al. 
2008); and (ii) high-efficiency coronagraphs such as the 
apodized pupil Lyot coronagraph (Soummer 2005 and ref- 
erences therein) or the achromatic 4-quadrant phase mask 
(Rouan et al. 2000; Mawet et al. 2006) to obtain optimal 
star extinction. GPI (Gemini Planet Imager) for Gem- 
ini South (Macintosh et al. 2006) and SPHERE (Spectro- 
Polarimetric High-contrast Exoplanet REsearch) for the 
ESO-VLT (Beuzit et al. 2006) are the two leading instru- 
ments of that category. These will both start operation in 
2011, along with HiGIAO (High-contrast Goronagraphic Im- 
ager for Adaptive Optics) for Subaru (Hodapp et al. 2008). 
They will aim at detecting exoplanets down to the Jupiter 
mass (Mjup) around nearby young stars by reaching con- 
trast values of 15 to 17.5 mag (10~® to 10~^) at angular 
separations of ~0.1". Both GPI and SPHERE will incor- 
porate diffraction limited integral field spectrographs (IPS) 
in the near infrared, allowing to obtain images simultane- 
ously at several wavelengths. SPHERE will also incorporate 
a differential imager named IRDIS (InfraRed Dual Imaging 
Spectrograph, Dohlen et al. 2008) that will provide simulta- 
neous images at two close wavelengths in either one of its 5 
different filter pairs over the Y to Ks bands. 

These instruments will allow to use different observ- 
ing strategies such as Spectral Differential Imaging (SDI, 
Racine et al. 1999) or Angular Differential Imaging (ADI, 
Marois et al. 2006) to obtain data which will be analyzed 
using advanced methods such as Spectral Deconvolution 
(Sparks & Ford 2002) or Localy Optimized Combination 
of Images (LOCI, Lafreniere et al. 2007) for IRDIS data. 
Specific methods of signal extraction have also been devel- 
oped within the SPHERE consortium (Mugnier et al. 2008; 
Smith et al. 2009) to be used in the data reduction pipeline 
of IRDIS. Data analysis methods are of extreme importance 
when it comes to the detection of faint objects in coron- 
agraphic images dominated by speckles. In particular, the 
precise estimation of the object fiux after applying these 
methods is critical for calibration of planetary mass ob- 
jects model atmospheres (Allard et al. 2001, 2003, 2007, and 
in preparation; Burrows et al. 2006; Ackerman & Mar ley 
2001; Tsuji 2005) based on effective temperature (T^eff) 
and surface gravity (log g), and the corresponding evolu- 
tionary models (Burrows et al. 1997; Chabrier et al. 2000; 
Baraffe et al. 2003; Fortney et al. 2005). The only currently 
available method for mass estimation, when no dynamical 



mass estimations are known for the object, consists in in- 
corporating photometric measurements and their error as- 
sessment into these physical models. Although spectroscopy 
is the method of choice for characterization, it may not be 
possible to obtain high-quality spectra of extremely faint 
sources (Vigan et al. 2008; Janson et al. 2010). 

In this paper we investigate the photometric limita- 
tions in high-contrast data obtained at different wavelengths 
with a dual-band imager like IRDIS, and we study how it 
translates in terms of characterization of planetary mass ob- 
jects. After reminding in Sect. 2 the origin of the speckle 
noise which is the fundamental limitation in high-contrast 
imaging and the methods to overcome it, we present in 
Sect. 3 the end-to-end simulations of IRDIS which were per- 
formed to obtain a realistic observing sequence. Section 4 
briefly describes the detection limits obtained with ADI and 
SDI-f ADI data analysis methods before studying the perfor- 
mances in terms of photometric accuracy. Finally, in Sect. 5 
we advocate a filter pair procedure for IRDIS and present 
an analysis of the characterizations which will be possible 
from aperture photometry. 



2 LIMITATIONS IN HIGH-CONTRAST 
IMAGING 

Detecting very faint planetary objects requires to obtain 
diffraction limited images with high-order AO systems in 
order to overcome the large contrast ratio between the star 
and the planet with coronagraphy. In high-Strehl ratio coro- 
nagraphic images, the factor that limits the accessible dy- 
namical range is the speckle noise (Soummer et al. 2007) 
induced by atmospheric phase residuals and instrumental 
quasi-static aberrations not corrected by the AO system. 
The quasi-static speckles are caused by the instrumental 
aberrations that slowly change during a long exposure. Tele- 
scope orientation, temperature variations or rotating optical 
elements cause small mechanical variations in the optical el- 
ements which make the speckle pattern evolve. Long time 
exposures are typically decomposed in a series of short expo- 
sure images of a few seconds during which the atmospheric 
residuals are averaged out, forming a smooth halo over which 
the quasi-static speckles are superimposed, because their co- 
herence time is much longer than the atmospheric residuals 
(Langlois et al. 1998; Macintosh et al. 2005; Hinkley et al. 
2007). To optimize AO performances and speckle rejection, 
observations can be performed in pupil stabilized mode, 
leading to a very high stability of the star Point Spread 
Function (PSF) and a slow rotation of the field of view dur- 
ing the observations, at a rate which depends on the star 
position in the sky. 

The speckle noise can be reduced by subtracting a ref- 
erence PSF from each science frame in order to remove the 
star halo and speckles, and possibly reveal a faint planetary 
object. This reference PSF can be obtained by observing a 
reference star taken in the same observing conditions (par- 
allactic angle and atmospheric conditions if possible) as the 
original target to reproduce a similar pattern of quasi-static 
speckles. This is a very time consuming task because the 
time spent on the reference star is equal to the one spent 
on the target to precisely match both PSF, and the aberra- 
tions between the two stars cannot be exactly reproduced. 
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The reference can also be built from the science frames by 
using either spectral or angular information. 

The SDI method has first been proposed by 
Racine et al. (1999) for faint companion detection. It 
has been extensively studied (Marois et al. 2000), and 
subsequently tested on sky with TRIDENT on the 
CFHT (Marois et al. 2005) and with NACO on the VLT 
(Lenzen et al. 2004). The technique relies on the fact that 
planetary objects have large intrinsic molecular features in 
their spectrum, while the host star has a relatively flat spec- 
trum. By taking simultaneously two images of a system at 
two close wavelengths located around one of these sharp 
features and subtracting them, the star contribution can be 
partially eliminated, and the planet signal revealed. SDI is 
most effective when used for detecting cool companions that 
show deep molecular absorption bands caused by H2O, CH4 
and NH3 at low Teff according to state of the art planetary 
mass objects atmosphere models. With carefully selected fil- 
ter pairs, a contrast of several magnitudes on the planet flux 
between the two filters can be obtained. However, the pres- 
ence of the molecular features expected for the detection of 
cool planetary companions should not be taken for granted 
since recent atmosphere models (Fortney et al. 2008), as 
well as observations of 2M1207b (Mohanty et al. 2007; 
Patience et al. 2010) and the IIR8799 planets (Marois et al. 
2008; Metchev et al. 2009; Janson et al. 2010), seem to show 
that non-equilibrium CO/CH4 chemistry could play an im- 
portant role in young low surface gravity objects. In partic- 
ular, the CII4 band head near 1.6 \ira could appear at much 
lower Te// than predicted by current atmosphere models. 
The SDI method is straightforward to implement: the im- 
ages taken at Ai need to be spatially rescaled to account 
for the spectral dependence of the PSF and subtracted from 
the images at Ao, with a possible amplitude correction factor 
to minimize the residual speckle noise. The main advantage 
of SDI is to significantly reduce the seeing halo, but it is 
intrinsically limited by speckle chromaticity and differential 
aberrations when going through two separate optical paths. 

The ADI method proposed by Marois et al. (2006) re- 
quires observations made in pupil-stabilized mode. It uses 
the field rotation to build an optimized reference PSF that 
contains very little signal from the planet. For each image 
li, a reference PSF is calculated using images taken before 
and images taken after, and for which a field rotation of at 
least 1.5 \/D has occurred in between. These images are 
then combined to eliminate the planet signal and produce a 
reference PSF that is subtracted from the image U. These 
operations are performed for all images in annuli of increas- 
ing radius. A thorough description of the complete procedure 
can be found in Sect. 5.2 of Marois et al. (2006). This tech- 
nique is essentially limited by the temporal evolution of the 
speckles which cannot be controled. The global efficiency of 
the ADI method is controled by the rotation rate of the field 
of view, which depends on the star declination, and by the 
angular separation, which constrains the actual motion on 
the detector. At the latitude of ESO-Paranal observatory 
(-24° 3' 38") and for a star at declination 5 = -45°, the 
field rotation varies between 0.006 deg s~^ at an hour angle 
of ±2h and 0.011 deg s~^ at an hour angle of Oh. This de- 
fines a strong constraint on the telescope time necessary to 
calibrate the speckles. 

Finally, the SDI and ADI methods can be efficiently 



Table 1. List of IRDIS filter pairs. 



Pair name 


Filter 


Filter 1 




Ao 




Al 










{\im) 




Y2Y3 


1.020 


20 


1.073 


20 


J2J3 


1.190 


25 


1.270 


25 


H2H3 


1.587 


30 


1.667 


30 


H3H4 


1.667 


30 


1.731 


30 


K1K2 


2.100 


20 


2.244 


20 



Table 2. Atmosphere models included in our library. 



Model 




log 


9 




(K) 


(dex) 


AMES-COND'' 


350 - 1300 


2.5 - 


6.0 


BT-SETTL'' 


1100 - 2300 


4.5 - 


5.5 


AMES-DUSTV^ 


1700 - 2500 


2.5 - 


6.0 


BSH'i 


700 - 2000 


4.5 - 


5.5 



=1 AUard et al. (2003) 
^ Allard et al. (2007) 
" Allard ct al. (2001) 
Burrows et al. (2006) 

combined to further reduce the speckle noise. SDI is first per- 
formed on short exposure images acquired simultaneously to 
remove the fast varying atmospheric residuals that have not 
been averaged out. ADI is subsequently applied on this set of 
data to combine the images with different angular positions 
of the field of view. 



3 END-TO-END SIMULATIONS 

A complete end-to-end model of SPHERE has been devel- 
oped to test the instrument performances and different data 
analysis methods. This model is a diffractive code written in 
IDL (Interactive Data Language) based on the CAOS (Code 
for Adaptive Optics System) problem solving environment 
(Carbillet et al. 2004) with a specific package developed for 
the SPHERE project (Carbillet et al. 2008). 

Realistic data cubes have been simulated to represent 
typical 4 hours exposure with IRDIS and an apodized Lyot 
coronagraph at different wavelengths where the star goes 
from -2 to ~f 2 hour angle. Every data cube is composed of 
144 images representing a cumulative 100 seconds esposure 
each, and several parameters are modified in the course of 
the simulation to take into account the variations of optical 
aberrations on a long timescale. The seeing and wind speed 
have been varied on ranges typical for the ESO-Paranal ob- 
servatory, respectively 0.85 ± 0.15" and 14.2 ± 4.6 m s~^. 
For each individual images, the AO-corrected atmosphere 
was simulated by a set of 100 decorrelated phase screens. 
The typical millisecond timescale of the uncorrected atmo- 
spheric residuals is not considered here: we assume that on 
a 100-second timescale, these residuals are averaged out and 
produce a smooth halo; only the correlated residuals with 
timescales longer than a few hundred seconds will remain. 
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Figure 1. l-cr noise levels after applying ADI and SDI+ADI data analysis methods in H band for high mass (left, FO at 10 pc, V = 2.7) 
and low mass (right, MO at 10 pc, V = 8.8) stars during a 4 h exposure time. The hatched area below 0.1" is covered by the opaque 
coronagraphic mask. The coronagraphic profile is calculated by the average of the coronagraphic image in annuli of increasing radius, 
and the different noise levels by the standard deviation of the residuals in the same annuli. All curves are normalized to the maximum 
of the PSF without coronagraph. 



From the instrumental point of view, variations of the beam 
shift as well as rotation of the entrance window, atmospheric 
dispersion corrector (ADC) and derotator have been trans- 
lated into wavefront error. Chromatic shifts associated to 
the ADC have been calculated from its optical design. Slow 
achromatic drifts, such as defocus and tilt, associated with 
temperature changes have been added. Finally, differential 
aberrations between the two filters in the differential imager 
have been taken into account: considering the prototypes 
of IRDIS DBI filters (Dohlen et al. 2008), 7.55 nm RMS of 
differential aberrations have been introduced. 

In these end-to-end simulations, the Fresnel propaga- 
tion of the wavefront is not considered. However, the overall 
impact of Fresnel propagation has been evaluated in sepa- 
rate simulations (not detailed here) where pre-coronagraphic 
and post-coronagraphic propagation effects have been sim- 
ulated. The main result is that while the region beyond AO 
cut-off (0.8" in If band) is mostly dominated by the pre- 
coronagraphic propagation effects, resulting in a loss of up 
to a factor 2 in contrast, the inner region is affected by a 
mix of both effects, resulting in a loss of at most 1.5. 

Four complete data cubes have been simulated, corre- 
sponding to the filter pairs Y2Y3, J2J3, H2H3 and K1K2 
of IRDIS (Table 1). The output of the diffractive code is a 
series of normalized coronagraphic and non-coronagraphic 
images of the star at the two wavelengths of the considered 
filter pair. A second code was used to create data cubes rep- 
resenting realistic planetary systems. For each star, 3 series 
of planets separated by 120° have been simulated at 0.2", 
0.5", 1.0", 1.5" and 2.0", taking into account the slow field 
rotation which is a function of the star elevation. The star 
was chosen at a declination of -45° and an hour angle of -2 h 
at the beginning of the simulated observation, representing 
a total field rotation of ~120°. 

To calculate realistic photometry we used standard Ku- 
rucz models (Kurucz 1979; Castelli & Kurucz 2003) for stars 
with spectral types regularly distributed from FO to MO at 
a distance of 10 pc (V = 2.7 to 8.8). For the planets we con- 
stituted a library of ~220 synthetic spectra including the 



AMES-Dusty models of AUard et al. (2001), the BT-Settl 
models of AUard et al. (2007), the AMES-Cond models of 
AUard et al. (2003) and the models of Burrows et al. (2006) 
with effective temperature ranging from Te// = 350 K to 
Te// = 2500 K, and surface gravity ranging from log g = 2.5 
to log g = 6.0. The steps in the grids of models are of 100 K 
in Te// and 0.5 in log g. We assume these models are comple- 
mentary, and Table 2 gives a list of the models with the Te// 
and log g ranges over which they are considered. For each 
filter pair we generated 66 data cubes with different com- 
bination of star and planet models to cover contrast values 
from 5 mag to 16.5 mag (--2 x 10"'' to ~10"^). 

In each data cube, sky contribution has been added to 
match typical values for ESO-Paranal observatory. Thermal 
background from the instrument was also included: the value 
is low in I, J and H bands (< 2 photon sec^^ pixeP^), 
while it becomes significant in K band (60 to 220 pho- 
ton sec~^ pixel"^). The code also accounts for the global 
throughput of the instrument and the atmospheric trans- 
mission, but does not consider OH lines variability. Finally, 
a realistic amount of noise for the IRDIS detector was in- 
cluded in the images: photon noise. Hat field noise (0.1%) 
and read-out noise (10 e^/read). 

The final output of the photometric code represents a 
4 h observation with IRDIS after standard cosmetic cor- 
rection and calibrations (dark, sky background and ther- 
mal background subtraction, flat field division, bad pixels 
correction). Due to the large number of parameters taken 
into account and the important computing time required 
for the simulation, only one data set representing a stan- 
dard case has been produced. This means that all generated 
data cubes present the same speckle pattern, and only differ 
by the photometric and noise values. 

ADI and SDI-I-ADI data analysis methods were then 
applied on all our simulated data cubes to attenuate the 
speckle noise. The ADI data analysis was implemented 
in IDL following the algorithm described by Marois et al. 
(2006): frames separated by 2.0 X/D were selected and com- 
bined in 5 annuli covering our simulated planets to produce 
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2 final images for data taken at Ao and at Ai. Tlie SDI data 
analysis was implemented in IDL using a custom routine 
for the precise spatial rescaling (L. Mugnier, private com- 
munication) based on zero-padding in both real and Fourier 
spaces. ADI was then applied on the subtraction of data 
taken at Ao and at Ai to produce the final SDI+ADI image. 



4 PHOTOMETRIC ACCURACY 

4.1 Noise level with ADI and SDI+ADI 

1-cr noise level were evaluated from the data products of 
the ADI and SDI-I-ADI data analysis methods by measur- 
ing the standard deviation of the residual speckle noise 
in annuli of increasing radius, normalized to the max- 
imum of the PSF without coronagraph. Although it is 
known from Goodman (1968), Aime & Soummer (2004), 
Fitzgerald & Graham (2006) and Smith et al. (2009) that 
the speckle noise statistics in AO-corrected images with and 
without coronagraph is not Gaussian. Marois et al. (2008) 
have shown that residual noise after applying ADI on ~20 
or more images is quasi-Gaussian. ADI was applied on our 
data cubes with more than 20 images in every case, so we 
will consider the residual noise as Gaussian and we will use 
standard 5-(t level for detection limits 

Fig. 1 shows the l-cr noise level for H2H3 filter pair 
in two different regimes: a high flux case corresponding to 
a high mass star (FO star at 10 pc, V = 2.7) and a low 
flux case corresponding to a lower mass star (MO star at 
10 pc, V = 8.8). The ADI noise level in flher H3 is not 
shown because it is at the same level as in filter H2. In 
high flux, the speckle noise attenuation is almost constant 
from 0.1" (inner- working angle of the coronagraph) to 2.5", 
both with ADI and SDI-I-ADI, allowing to reach a contrast 
of ~16 mag at 0.2" and more than 20 mag at 2.5". In low 
flux, the level of the background noise (sky and instrumental 
thermal emission, read-out noise) becomes limiting, and the 
attenuation reaches an almost constant level in SDI-f ADI 
at 1.5". The noise levels are similar in Y, J and H bands, 
but in K band, where the thermal emission is 10 to 15 times 
higher than in J or H band, the accessible contrast in low 
flux regime is limited at 15.5 mag. 

If we compare these 1-cr noise level to the ones derived 
for IRDIS Long Slit Spectroscopy (LSS) mode (Vigan et al. 
2008), we see that they are comparable between ADI and 
LSS. However, SDI+ADI clearly brings an improvement of 
2 to 3 magnitudes compared to LSS, justifying the need 
to quantify characterization capabilities in DBI mode for 
planets that are not detectable with LSS. 



4.2 Planet flux estimation 

We have estimated the signal of all planets detected at more 
than 5-cr with aperture photometry in a 2.44 \/D radius 
aperture. The aperture is calculated to take into account the 
effect of using a round aperture on square pixels. The mea- 
sured signal has been converted to a flux in phot s^^ m^^ 
with the formula: 



where Si is the measured signal on image i, Es is the encir- 
cled energy in the aperture, StcI is the telescope collecting 
surface, Tr is the transmission of the atmosphere, telescope 
and instrument, and ti is the exposure time for image i (100 s 
in our case). 

We consider that we are in a case where we know the 
value of the encircled energy Ee in the aperture. This value 
varies mostly with seeing conditions because the AO cor- 
rection will concentrate more energy in the PSF core when 
seeing improves. Moreover, we consider that the planet po- 
sition is known exactly to center the aperture on the planet 
PSF and avoid photometric error bias induced by inaccu- 
rate centering. Finally we also take into account the error 
induced by the field rotation. When the field rotates, it will 
slowly smear the planet PSF, especially at large angular sep- 
arations. The effect in our case is significant at separations 
larger than 1.0" because we simulated long exposures for 
individual images (100 s). In practice exposures for individ- 
ual images will typically last a few seconds to avoid detector 
saturation, reducing the effect of PSF smearing up to a few 
arcseconds. 

4.3 Photometric accuracy in ADI 

For our simulated test case, the planet flux has been evalu- 
ated in all data cubes after using ADI data analysis method. 
In each filter pair and for each simulated planet, 2 indepen- 
dent values are obtained at Ao and Ai . They are compared to 
the flux value introduced at the beginning of the simulation 
to evaluate the photometric error. Figure 2, left, illustrates 
the photometric performance as a function of wavelength 
and angular separation. The contours indicate the contrast 
value between the star and planet below which the photo- 
metric precision is better than 0.2 mag. Such a precision is 
necessary in order to be able to disentangle between differ- 
ent planet masses when comparing actual measurements to 
evolution nary models. We see two major effects: (i) the pho- 
tometric performance clearly depends on wavelength, and 
(ii) there are two different regimes depending on the posi- 
tion compared to the AO control radius. The first effect is 
directly related to the chromaticity of the PSF: in speckle- 
limited regime the noise attenuation is almost constant with 
angular separation compared to the coronagraphic profile, 
and the level of the coronagraphic profile linearly depends 
on wavelength. The second effect is related to the AO cor- 
rection inside the control radius. Inside that region we see a 
stabilization of the performance: 0.2 mag photometric pre- 
cision can be reached up to contrast of 10 to 11 mag (10~* 
to 4 X 10"*) from 0.2" to the AO control radius, which ex- 
tends from 0.5" in Y band to 1.0" in K band. Outside of the 
AO control radius, the photometric performance increases 
almost linearly with angular separation at all wavelengths 
to reach contrast values of 14 to 15 mag (2.5 X 10"" to 10"®) 
around 2.0". These numbers are given in the context of our 
simulated test case, but the general effects should be simi- 
lar for any data obtained with high contrast coronagraphic 
imagers. 

4.4 Photometric accuracy in SDI+ADI 



(1) 



Similarly to the noise level, using the SDI+ADI data anal- 
ysis method improves the photometric accuracy. However, 
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ADI / error 0.2 mag SDI+ADI / error 0.2 mag 



Y2Y3 J2 J3 H2 H3 K1 K2 Y2Y3 J2J3 H2H3 K1K2 




1.2 1.4 1.6 1.8 2.0 2.2 1.2 1.4 1.6 1.8 2.0 

Wavelength (/im) Wavelength (/.im) 



Figure 2. Magnitude difference between the star and the planet for which the photometric precision is better than 0.2 mag, as a function 
of wavelength and angular separation using ADI (left) and SDI+ADI (right) data analysis methods. The oblique dashod-line shows the 
AO control radius limit. 



using SDI+ADI will only provide an estimation of the dif- 
ferential flux of the planet between the 2 filters, contrary to 
ADI which provides an absolute measurement. To preserve 
the planet differential flux, the amplitude correction factor 
usually applied for SDI in the subtraction is taken equal to 
a fixed value of 1. The photometric error estimated with 
SDI+ADI follows the same variations as for ADI, but at 
higher contrast values. Figure 2, right, illustrates the photo- 
metric performance as a function of wavelength and angu- 
lar separation in SDI+ADI. The trends are similar to ADI 
alone, but the chromatic effect is less significant because 
the PSF chromaticity has been mitigated by the SDI part 
of the analysis. Compared to ADI alone, the contrast val- 
ues at which a 0.2 mag photometric error is reached are 1.5 
to 2.5 mag higher. At shorter wavelengths, in Y2Y3 filters, 
performances at separations larger than 1.0" decrease. This 
effect is related to the size of the aperture for photometry 
which is very small in Y band (4 pixels in diameter), and 
to the field rotation which has a strong effect on encircled 
energy at separations larger than 1.0" in Y band. Consider- 
ing shorter exposures for individual images where the field 
rotation is negligible would decrease the photometric errors 
in that particular case. 



4.5 Empirical photometric accuracy 

We hereafter combine the photometric accuracy obtained 
in ADI and SDI+ADI to define empirical photometric er- 
ror curves for each filter pair as a function of contrast. The 
photometric error curves as a function of contrast at each 
angular separation have been fitted with the empirically de- 
fined function: 

photorr = +P-i (2) 

where photcrr is the photometric error, c the contrast and 
(pi,P2,P3) the fitted parameters. This function approaches 
the measured points with a precision of ~1%. The fitting 
has been performed for ADI and SDI+ADI. To take into 
account the scattering of the error with the planet position 



Table 3. Contrast limit over which the differential photometric 
error in SDI+ADI becomes smaller than the photometric error in 
ADI. 



Separation Filter pair 

Y2Y3 J2J3 H2H3 K1K2 
(arcsec) (mag) (mag) (mag) (mag) 



0.2" 


6.8 


6.8 


6.5 


6.8 


0.5" 


6.0 


6.8 


8.3 


8.0 


1.0" 


10.8 


9.3 


8.0 


7.2 


1.5" 


12.8 


11.0 


11.0 


10.0 


2.0" 




11.5 


12.2 


12.2 



in the images, different cases have been considered at each 
separation, corresponding to the 3 different simulated planet 
positions: a standard case with an average photometric er- 
ror, an optimal case corresponding to the lowest estimation 
of the error and a pessimistic case corresponding to the up- 
per estimation of the error. These empirical photometric er- 
rors are plotted in Fig. 3 for the 4 simulated filter pairs. 
The amplitude of the error bars is defined by the optimal 
and pessimistic error curves described above. We assume 
that the photometric error in ADI is the same in the two 
filters of a pair, which is legitimate given the amplitude of 
the error bars. These empirical error curves lie in the same 
range as the expected photometric accuracy of other data 
analysis methods developed within the SPHERE consortium 
by Mugnier et al. (2008) and Smith et al. (2009). 

Table 3 gives for each filter pair and each angular sepa- 
ration the contrast value at which the photometric error in 
ADI becomes lower than the differential photometric error 
in SDI+ADI. These values give the contrast at which it be- 
comes more interesting in terms of photometric error to ob- 
tain a differential flux estimation. As explained in Sect. 4.4, 
aperture photometry in Y band is extremely sensitive to er- 
rors introduced by the position of the aperture or the field 
rotation because the aperture is very small. This is why in 
Y2Y3 pair at 2.0" there is no contrast limit between ADI 
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Figure 3. Empirical photometric error as a function of contrast in IRDIS filter pairs using either ADI or SDI+ADI data analysis 
methods. Errors bars have been represented only for a small set of data points. Their amplitude is defined by the optimal and pessimistic 
error curves described in the text. 



and SDI+ADI: for that particular case the flux estimation 
error is slightly better in ADI than SDI+ADI. 



5 PHOTOMETRIC CHARACTERIZATION 

In this section we evaluate characterization capabilities of 
IRDIS in imaging mode, i.e. how well the physical param- 
eters Te// and log g of the planets can be estimated from 
photometric measurements in different spectral bands. 

5.1 Characterization simulation 

To estimate the characterization capabilities of IRDIS, we 
performed a new simulation using as input the 5-cr detec- 
tion limits obtained from Sect. 4.1 and the empirical error 
curves obtained in Sect. 4.5. The goal of the simulation was 
to test the efficiency of all filter pair sequences for char- 
acterization at different stellar magnitudes and for a large 
number of planetary atmosphere models. These simulations 
are based on current state of the art atmosphere models. Al- 
though these models will clearly evolve with new detections 
in the future, they allow to test the expected performances 
of IRDIS, as well as to estimate the intrinsic errors of our 



signal extraction and comparison to models. It has been per- 
formed for all stellar types and atmosphere models included 
in our library (see Table 2). 

For the simulation we assume that a same planetary 
system is observed with different filter pairs in a given or- 
der. For each possible combination of parameters (filter pairs 
sequence; star magnitude; angular separation; planet atmo- 
sphere model) we proceed as follow: the star and planet 
fluxes are calculated in the filters of the first pair; if the 
planet is not detectable (considering the 5-(t detection limit), 
simulation for that combination of parameters is stopped; if 
it is indeed detectable, a photometric measurement is ob- 
tained. Depending on the contrast between the planet and 
the star in each filter, different informations are obtained: 2 
direct photometric measurements if the planet is detectable 
with ADI in both filters, a differential fiux measurement if 
the planet is only detectable with SDI+ADI or a direct and 
a differential measurement if the planet is only detectable 
in ADI in one of the filters. Once the flux measurements 
are obtained, the photometric error is determined from the 
empirical error curves and added to the measured values 
to obtain lower and upper limits to the planet flux. Mod- 
els that can correspond to these limits are then searched in 
our models library. If only one model corresponds, we stop 
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iterations considering that the planet has been fuUy charac- 
terized within the limits of the grid of atmosphere models 
that is used. If many models match those limits, we switch 
to the next filter pair in the sequence and the same process 
is started again. In a given sequence, each filter will bring 
some additional information that will help to find the ap- 
propriate atmosphere model and constrain the values of Te// 
and log g. At the end of the filter pairs sequence, 4 distinct 
outcomes are possible: 

• Non-Detection (ND): the planet is not detected in the 
first filter pair; the sequence is stopped. 

• Non-umque Characterization (NC): the planet is 
detected at least in the first filter pair of the sequence; 
at the end of the sequence many models match the flux 
measurements and they have different values of T^ff and 
log g- 

• Te// Characterization (TC): the planet is detected at 
least in the first filter pair of the sequence; at the end of the 
sequence many models match the flux measurements and 
they all share the same value for Te// but not for log g. 

• Full Characterization (FC): the planet is detected at 
least in the flrst filter pair of the sequence; at the end of the 
sequence only one model matches the flux measurements, 
which means that the Te// and log g values have been de- 
termined. 

The TC and FC are considered within the limits of the 
atmosphere models grid, which is 100 K in Te// and 0.5 
in log g, i.e. the full characterization corresponds to the 
determination of Te// and log g with an error equal to the 
limits of the models grid. Another important point is that 
in practice the photometric error will have to be estimated 
from the science data itself. Although, we consider here an 
ideal case where the photometric error is known, other data 
analysis methods such as the one proposed by Smith et al. 
(2009) will allow direct estimation of the error from the data, 
with a precision that it compatible with the results presented 
here. 



5.2 Filter pair sequence analysis 

The output of this simulation allows to determine the most 
significant filter pairs sequence for characterization, i.e. the 
sequence which maximizes the number of characterizations. 
The strategy is to progressively build an optimal sequence 
by adding each time the filter pair that increases the most 
the number of characterizations. All possible filter pair se- 
quences have been systematically tested to find the one that 
maximizes the number of Te// characterizations and full 
characterizations (TC+FC), as a function of stellar magni- 
tude. The simulation shows the following important results 
(see Table 4 for details): 

• H2H3 is the filter pair which minimizes the number 
of non-detections, reflecting the fact that CH4 absorption 
band near 1.6 \im in cool substellar objects spectra is the 
optimal feature for their detection. However it should be re- 
minded that this speciflc spectral feature might not always 
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Figure 4. Proportion of models in our library characterized by 
adding filter pairs from the optimal filter pair sequence for 4 stel- 
lar magnitudes. The error bars are given by the optimal and pes- 
simistic photometric error curves. 



be present in lower mass objects as we mentionned in Sect. 2. 

• When adding a new filter pair to a sequence, adding 
Y2Y3 or J2J3 will increase the number of characterized 
models more than H3H4 or K1K2. This result does not 
depend on the stellar magnitude. 

• The scattering of the flux error between the pessimistic 
optimal photometric error curves has no major influence on 
the number of characterized models when more than one 
filter pair is used. In particular, we see the error bars of the 
different fiux cases do not significantly overlap, confirming 
that the chosen sequence is appropriate for all cases of flux. 

This last result is particularly important as it advocates 
for a given sequence of filter pairs and allows to set general 
priorities on the filter pairs for characterization. When there 
are no a priori assumptions on the nature of the objects, the 
priorities are defined from highest to lowest as follow: 

0. H2H3 

1. Y2Y3 / J2J3 

2. H3H4 / K1K2 

The H2H3 - Y2Y3 - J2J3 ~ H3H4 - K1K2 wiU be referred to 
as the optimal sequence from now on, and we will only con- 
sider the standard empirical errors since the differences with 
the other error curves are small. Assuming the use of the 
optimal sequence. Fig. 4 represents the proportion of char- 
acterized models (TC+FC) from our library when new filter 
pairs are added. The trends are identical for the 4 different 
stellar magnitudes. When using only H2H3 the proportion 
of characterized models is comprised between 30 and 80%, 
and the error bars are of ~10% for bright stars. Adding a 
second filter pair greatly improves the proportion of charac- 
terized models which is above 70% for all stellar magnitudes. 
Adding more filter pairs confirms this trend and tends to re- 
duce the error bars to less than 5% for all magnitudes. The 
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Table 4. Analysis of the filter pairs sequences. 
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main conclusion is that most of the information for char- 
acterizing any given model is obtained using 2 filter pairs 
around low mass stars when the contrast is favorable, and 3 
filter pairs around high mass stars for which the contrast is 
more challenging. 

5.3 Lowest estimations of Te// 

In the previous section we set priorities for characterization 
on the different filter pairs of IRDIS. We will now detail the 
lowest values of Te// that IRDIS will be able to characterize 
as a function of stellar magnitude and angular separation. 
Figure 5 gives the smallest values of Te// which have been 
characterized when using 1 to 3 filter pairs from the optimal 
sequence. Colder planets can be detected, but we were not 
able to find the appropriate values of Te// and log g. When 
using only H2H3, planets with Te// down to 900 K should 
be characterized at an angular separation of 0.2" from high 
mass bright stars and 700 K from lower mass stars. Adding 
a second filter pair considerably improves these results by 
200 K, while adding a third pair confirms these limiting 
values. 

With the considered data analysis methods and accord- 



ing to the evolutionary models from Baraffe et al. (2003) for 
the COND atmosphere models, we can estimate that in a 
very young system of 10 Myr, we should be able to charac- 
terize a planet of 1 Mjup with H2II3 at separations larger 
than 0.5" around a low mass star (MO at 10 pc) where the 
star-planet contrast is favorable, but only further than 2.0" 
around a high mass star (FO at 10 pc) where the contrast 
difference is larger. With two filter pairs, the limit would 
be 0.2" around a low mass star and 1.0" around a high 
mass star. For older systems, only planets of a few masses 
of Jupiter could be characterized. At 100 Myr, a Jupiter 
mass planet would remain out of reach for characteriza- 
tion with H2H3 filters around a high mass star, and only 
at separations larger than 1.5" around a low mass star. At 
this age, the Te// limits of 700 K and 500 K which can be 
reached at small angular separation around high and low 
mass stars would respectively correspond to planets with 
masses of ~6.5 Mjup and ^3 Mjup. Using improved signal 
extraction methods providing more accurate photometry of 
the companion would certainly push down those limits. 
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Figure 6. Distribution of tlic errors on T^eff and log g for non-unique characterizations with H2H3 filter pair around a high mass star 
(FO at 10 pc, plain line) and a low mass star (MO at 10 pc, dashed line) cases. 



5.4 Study of the non-unique characterizations 

NCs are the cases where several models correspond to the 
flux measurements in all filter pairs with which they are 
detected. From these remaining models, it is possible to de- 
termine if a combination of (Te//; log 17) is more represented 
than others, making this combination the most probable val- 
ues of Te// and log g. If several combinations are counted an 
equal number of times, an average value and an error can be 
determined for the values of Te// and log g. In any case, the 
error is at least equal to the steps in the grids of models. The 
estimation of the most probable values for Te// and log g has 
been performed for all non-uniquely characterized models at 
all simulated angular separations and magnitudes. 

Figure 6 shows an histogram of the errors on Te// and 
log g when using H2H3 filter pair for high mass (FO at 10 pc) 
and low mass (MO at 10 pc) stars. NCs are mostly domi- 
nated by errors on the determination of log g. In particular 
we see that around a low mass star where the contrast is 
more favorable, the proportion of models with no error on 
Te// increases by '-^^20%, while the errors on log g keep the 
same distribution. Adding other filter pairs improves the de- 
termination of Te// for the non-unique charaterizations, in 
particular around high mass (brighter) stars, reaching more 
than 95% for all stellar magnitudes. The determination of 
log g is also improved, but even when using 3 filter pairs 
the number of cases where the error is less or equal to 0.5 
never reaches more than 85% around a high mass star or 
95% around a low mass star. 



5.5 Impact of errors on Te// and log g 

The influence of errors on the determination of Te// and 
log g on the determination of the planet mass can be stud- 
ied using evolutionnary models such as those published by 
Baraffe et al. (2003) for the COND atmosphere models. The 
preliminary version of the SPHERE target list (S. Desidera, 



private communication) was used to deflne a standard young 
test case based on age considerations. The average age for 
targets younger than 100 Myr is 44 ± 20 Myr and the av- 
erage error on the determation of the target age is 30 Myr. 
Using these values we can define two test cases considering a 
planet of 2 Mjup aged of 44±30 Myr orbiting at 5 A.U. from 
MO and FO stars at 10 pc. According to the evolutionary 
models from Baraffe et al. (2003), such planets should have 
Te// = 516 K and log g = 3.54 dex, resulting in a contrast of 
11.9 mag and 15.6 mag in H band respectively around the 
MO and FO stars. Considering the results from Sect. 5.3 and 
5.4, the planet around the FO star cannot be characterized 
with IRDIS using one filter pair, while the planet around 
the MO star is close to the measured limit. 

The expected spectra of these planets were introduced 
in our simulation to test the accuracy of extracting flux in- 
formation and inversely deriving physical parameters, and 
to estimate their mass using evolutionary models. The re- 
sults are presented in Fig. 7: the areas covered by the values 
of Te// and log g are shown as rectangular boxes on the pre- 
dicted isochrones for both planetary systems. The bin size 
of the atmosphere model grid is also plotted for reference. 
For each case, the planet mass is estimated by selecting all 
the isochrones of the masses that cross the error box and 
weighting them by the integral of the isochrone inside the 
box. The isochrone that has the largest intersection with the 
error box is supposed to be the most likely (or the average 
if several isochrones have the same integral) . The upper and 
lower limits of the mass are given by the highest and lowest 
mass isochrones that cross the box. 

Around a low mass star, the parameters Te// and log g 
are estimated with an accuracy close to the one given by the 
atmosphere model grid Mjup), leading to an estima- 

tion of 1.9lj Q Mjup. Around a high mass star, the planet is 
very close to the detection limit at 0.5", resulting in a poor 
estimation of both Te// and log g: the important photomet- 
ric error in H2H3 leads to a very large uncertainty on log g 



Exoplanet characterization using ADI and SDI 11 



H2H3 



2.0 Mj @ 0.5" / 40±30 Myr 




0.5 1.0 1.5 

Angulare separation (as) 

H2H3+Y2Y3 




0.5 1.0 1.5 

Angulare separation (as) 

H2H3+Y2Y3+J2J3 



0.5 1.0 1.5 

Angulare separation (as) 



2.0 




2.0 



Figure 5. Smallest value of Teff which can be characterized as 
a function of star magnitude (y-axis) and angular separation (x- 
axis) when using 1, 2 or 3 filters from the optimal sequence. 



(4.33 ± 1.23). The mass of the planet is then estimated to 
l.llg'g Mjup. In this case, the planet mass is stronly un- 
derestimated, and the large uncertainties on T^ff and log g 
lead to a large upper limit for the mass estimation. The ori- 
gin of the large offset on the determination of log g is still 
uncertain, and further simulations are still needed to investi- 
gate thoroughly the complete parameter space. In particular 
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Figure 7. Isochrones for the COND planetary atmosphere mod- 
els covering an age of 40 ± 30 Myr used for the determination of 
the mass of hypothetical 2 Mjup planets orbiting at 5 A.U. from 
MO and FO stars at 10 pc. The error boxes defined by the possible 
values for Teff ^^'^ log g of both planets are respectively repre- 
sented by dotted and dashed rectangles. The planet mass derived 
from the models and error box is displayed at the bottom left 
corner of each error box (see text for an explanation on how the 
planet masses are derived). The position of the planet predicted 
by the evolutionary models is represented by the star symbol, and 
the error box defined by the atmosphere models grid precision is 
given for reference around that position by a plain rectangle. 



we see in this case that the age uncertainty unquestionably 
increases the uncertainty on the mass estimation by increas- 
ing the number of isochrones crossing the possible values of 
Te// and log g. Such simulations would greatly benefit from 
updated homogeneous atmosphere models grids covering a 
large span of Te// , log g and age. 



6 CONCLUSIONS 

Next generation instruments for ground based exoplanet di- 
rect imaging such as SPHERE and GPI will provide data 
intrinsically limited by the speckle noise. This noise needs 
to be attenuated using a posteriori data analysis methods, 
such as spectral and angular differential imaging. In this pa- 
per, we have quantified the exoplanet characterization capa- 
bilities of IRDIS, the differential imager of SPHERE, using 
photometric and differential photometric information. 

The photometric performances have been evaluated 
with aperture photometry on the detectable planets as a 
function of contrast and wavelength for a standard test case. 
In particular we have shown that the photometric perfor- 
mance strongly depends on wavelength because of the PSF 
chromaticity and on the position with respect to the AO con- 
trol radius. With ADI, a photometric accuracy of 0.2 mag is 
reached inside the AO control radius for contrast values of 
10 to 11 mag between the star and planet, while at larger ra- 
dius the precision can be reached for contrasts up to 15 mag. 
With SDI+ADI, the overall photometric performance is im- 
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proved, increasing by 1.5 to 2.5 mag the contrast at which 
a 0.2 mag precision is reached. 

We have defined empirical photometric accuracies for 
IRDIS in its different filter pairs, which were used to test 
the characterization capabilities of IRDIS with all possible 
combinations of filter pairs. Priorities for characterization 
have been set on the different filter pairs by finding the 
pairs that maximize the number of possible characteriza- 
tions in various conditions. We showed that when there is 
no a prion knowledge on the planet, the filter pairs Y2Y3 
and J2J3 allow a larger number of characterizations than 
the pairs H3H4 and K1K2. Then we showed that using filter 
pair II2II3, it will be possible to characterize planets with 
Te// — 900 K around high mass stars at small angular sep- 
aration and Te// — 700 K around lower mass stars. Adding 
Y2Y3 and J2J3 filter pairs allows to decrease the charac- 
terizable Te// by 200 K at all separations and for all stellar 
magnitudes, while considerably decreasing the number of 
non-characterizations for warmer planets. 

Finally, we showed that non-unique characterizations, 
i.e. planets for which the Te// and log g values could not be 
determined exactly (within the limits of the grid of atmo- 
sphere models) are mostly dominated by errors on the de- 
termination of log g. In particular, we showed that around 
low mass stars, where the contrast is more favorable, the 
determination of Te// is largely improved, while the errors 
on log g remain identical around a high mass star. Consid- 
ering evolutionary models, and including typical ages from 
the future SPHERE target list, we showed that such errors 
on the determination of a low mass planet (2 Mjup) results 
in a large uncertainty around a high mass bright star, but 
is very close to the limits fixed by the models grid around a 
fainter low mass star. 

With this work we showed that IRDIS, the dual-band 
imager of SPHERE, should be able to fulfill the goal set 
for a high-contrast imager, that is the ability to detect and 
characterize planetary companions down to the Jupiter mass 
around nearby young stars. Similar developments could also 
be performed for an Integral Field Spectrograph in the 
future, allowing to quantify precisely the performances of 
SPHERE in the near-infrared, and to work on the aspect of 
characterization strategy for the detected objects. 
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